Transient Solutions of Markov Processes by Krylov Subspaces
نویسندگان
چکیده
In this note we exploit the knowledge embodied in innnitesi-mal generators of Markov processes to compute eeciently and economically the transient solution of continuous time Markov processes. We consider the Krylov subspace approximation method which has been analysed by Y. Saad for solving linear diierential equations. We place special emphasis on error bounds and stepsize control. We discuss the computation of the exponential of the Hessenberg matrix involved in the approximation and an economic evaluation of the Pad e method is presented. We illustrate the usefulness of the approach by providing some application examples. Approximation des solutions transitoires d'un processus de Markov dans des sous-espaces de Krylov R esum e : Nous nous int eressons dans cet article au calcul du r egime tran-sitoire d'un processus de Markov. Un tel probl eme n ecessite le calcul d'une exponentielle de matrice. En prenant en consid eration certaines propri et es inh erentes aux mod eles markoviens, nous examinons l'approximation dans un espace de Krylov propos ee par Y. Saad. Cela nous permet de pr eciser certaines caract eristiques th eoriques et de d eenir un algorithme able et ef-cace. Nous pr esentons des r esultats num eriques attestant de la validit e de la m ethode. Mots-cl e : Cha^ nes de Markov, exponentielle d'une matrice, espaces de Krylov, m ethode d'Arnoldi.
منابع مشابه
2nd International Workshop on the Numerical Solution of Markov Chains 1 Preconditioned Krylov Subspace Methods for the Numerical Solution of Markov Chains
In a general projection technique the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of techniques based on this principle is that of Krylov subspace methods which utilize subspaces of the form spanfv; Av; ::::;A m?1 vg. This general principle can be used to solve linear systems and eigenvalue problems whic...
متن کاملModel Checking Markov Chains Using Krylov Subspace Methods: An Experience Report
The predominant technique for computing the transient distribution of a Continuous Time Markov Chain (CTMC) exploits uniformization, which is known to be stable and efficient for non-stiff to mildly-stiff CTMCs. On stiff CTMCs however, uniformization suffers from severe performance degradation. In this paper, we report on our observations and analysis of an alternative technique using Krylov su...
متن کاملComputing Approximate (Symmetric Block) Rational Krylov Subspaces without Explicit Inversion
It has been shown that approximate extended Krylov subspaces can be computed –under certain assumptions– without any explicit inversion or system solves. Instead the necessary products A−1v are obtained in an implicit way retrieved from an enlarged Krylov subspace. In this paper this approach is generalized to rational Krylov subspaces, which contain besides poles at infinite and zero also fini...
متن کاملKrylov Subspace Methods in Linear Model Order Reduction: Introduction and Invariance Properties
In recent years, Krylov subspace methods have become popular tools for computing reduced order models of high order linear time invariant systems. The reduction can be done by applying a projection from high order to lower order space using bases of some subspaces called input and output Krylov subspaces. One aim of this paper is describing the invariancies of reduced order models using these m...
متن کاملLocally Optimal and Heavy Ball GMRES Methods
The Generalized Minimal Residual method (GMRES) seeks optimal approximate solutions of linear system Ax = b from Krylov subspaces by minimizing the residual norm ‖Ax − b‖2 over all x in the subspaces. Its main cost is computing and storing basis vectors of the subspaces. For difficult systems, Krylov subspaces of very high dimensions are necessary for obtaining approximate solutions with desire...
متن کامل